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Abstract 

We consider energy minimizing configurations of a nematic liquid crys¬ 
tal around a spherical colloid particle, in the context of the Landau- 
de Gennes model. The nematic is assumed to occupy the exterior of a 
ball B ro , and satisfy homeotropic weak anchoring at the surface of the 
colloid and approach a uniform uniaxial state as |x| —> oo. We study 
the minimizers in two different limiting regimes: for balls which are small 
ro L 2 compared to the characteristic length scale is, and for large 
balls, ro L 5 . The relationship between the radius and the anchor¬ 
ing strength W is also relevant. For small balls we obtain a limiting 
quadrupolar configuration, with a “Saturn ring” defect for relatively strong 
anchoring, corresponding to an exchange of eigenvalues of the Q-tensor. 
In the limit of very large balls we obtain an axisymmetric minimizer of 
the Oseen-Frank energy, and a dipole configuration with exactly one point 
defect is obtained. 


1 Introduction 

Liquid crystals are well-known for their many applications in optical devices. 
The rod-like molecules in a nematic liquid crystal tend to align in a common 
direction: the resulting orientational order produces an anisotropic fluid with 
remarkable optical features. This anisotropy also makes it highly interesting 
to use nematic liquid crystals in colloidal suspensions. Immersion of colloid 
particles into a nematic system disturbs the orientational order and creates 
topological defects, which enforce fascinating self-assembly phenomena [23, 19], 
with many potential applications [26, 22]. This sensitivity to inclusion of small 
foreign bodies also has promising biomedical applications [32]: for instance, new 
biological sensors could detect very quickly the presence of microbes, based on 
the induced change in nematic order [28, 13, 14]. 

In the present paper we investigate the structure of the nematic order around 
one spherical particle, with homeotropic [i.e. normal) anchoring at the particle 
surface, and uniform alignment far away from it. The homeotropic anchoring 
creates a topological charge. This charge must be balanced in order to match 
the uniform alignment at infinity, which is topologically trivial. Therefore one 
expects to observe singularities. 

This particular problem is a crucial step in understanding more complex 
situations, and it has received a lot of attention in the past two decades [31, 16, 
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17, 29, 24]. These works rely on heuristically supported approximations, and 
numerical computations. They point out two possible types of configurations, 
with “dipolar” or “quadrupolar” symmetry (related to their far-held behavior 
[30, § 4.1]). In a dipolar configuration the topological charge created by the 
particle is balanced by a point defect, while in a quadrupolar configuration it is 
balanced by a “Saturn ring” defect around the particle. 

The aforementioned works use either Oseen-Frank theory [31, 16, 17, 29] or 
Landau-de Gennes theory [24] to describe nematic alignment. In Oseen-Frank 
theory, the order parameter is a director held n(x) £ S 2 , which minimizes an 
elastic energy. One drawback of that model is that line defects have infinite 
energy. In particular, the energy of a quadrupolar configuration with Saturn 
ring defect has to be renormalized. Moreover, Oseen-Frank theory only accounts 
for uniaxial nematic states: it assumes local axial symmetry of the alignment 
around the average director. On the other hand Landau-de Gennes theory 
involves a tensorial order parameter that can also describe biaxial states, in 
which the local axial symmetry is broken. This is the model that we will be 
using here. 

The order parameter in Landau-de Gennes theory is the so-called Q-tensor, 
which belongs to the space 


So := {Q G M 3 (R): Qij = Q jit tr(Q) = 0} , (1) 


of symmetric traceless 3x3 matrices. The eigenvectors of Q represent the 
average directions of alignment of the molecules, and the associated eigenvalues 
measure the degree of alignment along these directions. Uniaxial states are 
described by Q-tensors with two equal eigenvalues, which can be put in the 
form 


Q = s 



s G R, n G S 2 . 


Biaxial states correspond to the generic case of a Q-tensor with three distinct 
eigenvalues. 

The configuration of a nematic material contained in a domain fl C R 3 is 
described by a map Q: Q —> So- At equilibrium it should minimize the free 
energy functional 


HQ) = 



|VQ| 2 + /(Q) 


dx + T S (Q). 


( 2 ) 


Here J~ S (Q) is a surface energy term which depends on the type of anchoring 
(see (6) below), and the bulk potential f(Q)> 0 is given by 


f{Q) = - | tr(Q 2 ) - jj tr(Q 3 ) + | tr(Q 2 ) 2 + C 0 , (3) 

for some material-dependent constants a > 0, b, c > 0. The constant Go is 
chosen to ensure that min / = 0. The set of Q G Sq minimizing the potential 
(3) obviously plays a crucial role. It consists exactly of those Q-tensors which 
are uniaxial, with fixed eigenvalues: 


74 := {/ = 0 } 



(4) 
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where s* = (b + \Jb 2 + 24ac)/4c > 0. 

We are interested here in the nematic configuration around a spherical par¬ 
ticle: 

n = a ro :=R 3 \B ro , 

where > 0 is the particle radius. We impose uniform W*-valued conditions at 
infinity 


lirn Q(x) = Qoo := s* ( e 2 ® e 2 - \l ) , e z = (0,0,1). (5) 

|x|—>oo \ 6 J 

At the particle surface, weak radial anchoring is enforced through the surface 
term J~ s in the free energy functional (2). This surface contribution is given by 

MQ) = Z f \Qs-Q\ 2 dA, ( 6 ) 

Z J dB r(J 

where W > 0 is the anchoring strength, and Q s is the W*-valued radial map 

Qs := s* (e r ® e r - = ]^f ( 7 ) 

Denoting by v the exterior normal to fl ro , the corresponding boundary condi¬ 
tions are 

Wl^v = ® s ~® for |ar| = r 0 . (8) 

We also include in this description the case of strong anchoring, corresponding 
to W — +oo and Dirichlet boundary conditions 

Q = Q s for \x\ = r 0 . (9) 


In every case, the Euler-Lagrange equations 

\Q \ 2 + c \Q\ 2 Q, (10) 

are satisfied in T>'(il ra ; iSq) by any equilibrium configuration. 

The existence of minimizers of the free energy functional (2) with the uniform 
far-field condition (5) can be obtained if we replace the pointwise condition (5) 
with the integrability condition 


LAQ = V/(Q) = —aQ -MQ 2 - 




dx < oo. 


( 11 ) 


The Euler-Lagrange equations (10) can then be used to see that the strong 
condition (5) is in fact also satisfied. After establishing this existence result in 
Section 2, we turn to studying the two asymptotics regimes of “small” or “large” 
colloid particle. 

According to the numerical computations in [29, 24], small particles favor 
quadrupolar configurations with a defect ring, while large particles favor dipolar 
configurations with a point defect. In the present paper we obtain rigorous 
justifications of these observations. In the small particle regime we also provide 
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exact information on the radius of the defect ring, for which the values computed 
in [31, 16, 17, 29, 24] did not agree. 

We determine the size of the particle based on the ratio r^/L, and consider 
the limits as the ratio tends to zero and to infinity. In either limiting case, the 
ratio rgW/L, which gauges the anchoring at the particle surface, will also enter 
into the description of the limiting configuration. 


The small particle limit. In Section 3 we investigate the small particle 
regime, and prove: 


Theorem 1. Consider (for any ?’ 0 , W or L) a map Q, finite-energy solution 
of (10)-(8)-(ll). Let w £ (0, Too] be the effective “limiting anchoring strength 
Then, as 



( 0 ,«>)> 


the rescaled maps x H > Q(tqx) converge to 


Qo s * 


w 1 

3 + wr 3 




+ s*(l 


w 1 

1 + w r 


) 




locally uniformly in Q. i. 

An interesting feature of Theorem 1 is the explicit form of the limit: it 
provides a very precise description of the quadrupolar configurations. The Sat¬ 
urn ring defect appears as a discontinuity in the principal eigenvector of Q(x), 
the Q-tensor passing through a uniaxial state as eigenvalue branches cross via 
an “eigenvalue exchange” mechanism [24]. The ratio of the ring radius to the 
particle size is thus found to be, for w > \/‘i, the solution r > 1 of 


r 


3 


w 

- - r 

1 + w 


w 

3 + w 


For the strong anchoring w = oo, its value is r « 1.47. As the anchoring 
strength decreases, the ring shrinks until it becomes a surface ring for w = y/3. 
At very weak anchoring w < y/3 there is no defect ring anymore. (See Figure 1 
in Section 3 for the ring location at various values of w.) This description is 
consistent with [31, 16, 17, 29, 24], with the significant improvement of providing 
exact values for the relevant quantities. 

Away from the ring defect, the limiting map Qo(x) is everywhere biaxial in 
fix: the small particle limit r^/L —> 0 does not correspond to the unit director 
Oseen-Frank model. Indeed, it is well known (see [25], and [8] for more general 
energy functionals) that minimizing S 2 -valued maps cannot have line defects. 


The large particle limit. The large particle regime is more delicate to anal¬ 
yse. We restrict ourselves to minimizers of the free energy, and to strong an¬ 
choring - that is, W = oo. For the rescaled maps x K > Q(rox), the regime 
rfj/L 1 corresponds to the vanishing elastic constant limit studied in [18, 20]. 
There, the authors prove convergence to a W*-valued map whose director is an 
S 2 -valued minimizing harmonic map. It is well-known that such a map has a 
discrete set of singularities [25], and that these defects carry topological degrees 
±1 [5]. In particular, the results of [25, 7] ensure that a “Saturn ring” defect 
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cannot be observed in the large particle limit, regardless of the anchoring con¬ 
dition on the particle surface, and thus very different behavior may be expected 
in the large particle regime than was observed for small particles. 

Since the strong radial anchoring imposes a degree +1 near the particle, 
while the uniform far-held condition imposes a zero degree at infinity, there 
must be at least one point defect of degree —1. However the number of defects 
of a minimizing map does not necessarily correspond to the minimal number of 
defects required by the topology [10]. In our case we expect that there is exactly 
one defect, as predicted by [17, 29, 24]. Since determining the exact number 
of defects is a very difficult question in general, we restrict ourselves to axially 
symmetric configurations: we impose invariance under any rotation of vertical 
axis, and that eg (horizontal unit vector orthogonal to the radial direction) be 
everywhere an eigenvector of the Q-tensor. This natural symmetry assumption 
seems to be supported by the numerical pictures in [24]. 

In the limit we will therefore obtain an axially symmetric S 2 -valued harmonic 
map. Such maps have been studied in [9, 12]. They are analytic away from a 
discrete set of defects on the 2 -axis. For very particular symmetric boundary 
data, it can be deduced from rearrangement inequalities that the number of 
defects matches the topological degree [9, Theorem 5.1]. This result does not 
apply to our case, but - using different arguments - we nevertheless manage to 
show that there is exactly one defect, thus justifying the dipole configuration 
predicted by [17, 29, 24], More precisely, in Section 4 we prove: 

Theorem 2. Let Q minimize the free energy (2) among axially symmetric 
maps satisfying the boundary conditions (9)-(ll). Then, as r^/L goes to +oo, 
a subsequence of the rescaled maps x <—> Q{r$x) converges to a map 


Q*(x) = s*(n(x) (g) n(x) — 1/3), 


locally uniformly in Hi \ {po}- Here n minimizes the Dirichlet energy in Hi, 
among axially symmetric S 2 -valued maps satisfying the boundary conditions 


n = e r on dB\, and 


K ) 2 


In i 


--P^dx< oo, 


and n is analytic away from exactly one point defect po, located on the axis 
of symmetry. 


The core of Theorem 2 is verifying that the minimizing harmonic map n 
admits at most one defect. We achieve this by investigating the topology of the 
sets {ii 3 > 0 } and {713 < 0 } where n points “more upward” or “more downward”. 
Using basic energy comparison arguments and the analyticity of minimizers 
away from the 2 -axis, we show that these sets are connected. Merging this 
with the observation that defects correspond to “jumps” between upward- and 
downward-pointing n, we conclude that there cannot be more than one defect. 

We note that there are very few cases in which the number of defects is 
actually known to match the topological degree. This is true in a ball with 
radial Dirichlet boundary conditions, because then the energy of the radial map 
can be explicitly computed and seen to coincide with a general lower bound [5], 
and this is true also for geometries close enough to the radial one [11], We do 
expect that our minimizers in the axially symmetric class are actually minimizers 
in the general (nonsymmetric) case, but this remains an open question. 


5 



The plan of the paper is as follows. In Section 2 we prove the existence of 
minimizers and some basic properties. In Section 3 we investigate the small par¬ 
ticle regime and the quadrupolar “Saturn ring” configurations. In Section 4 we 
study the large particle regime and the associated axially symmetric harmonic 
map problem. 

Acknowledgements: Part of this work was carried out while XL was visiting 
McMaster University with a “Programme Avenir Lyon Saint-Etienne” doctoral 
mobility scholarship. He thanks the Mathematics and Statistics Department of 
McMaster University for their hospitality, and his doctoral advisor P. Mironescu 
for his support and advice. 

Notations: We will use cylindrical coordinates (p, 9 , z) defined by 
X\ = pcos9 , X2 = psin0, X3 = z, 
and the associated orthonormal frame (e p ,eg,e z ), where 

e p = (cos 0, sin 0,0), eg = (— sin$, cos 9, 0). 

We will also use spherical coordinates (r, 9 , <p) with (r, ip) defined by 

p = r sin ip, z = r cos ip. 


2 Existence and first properties of minimizers 


We start by remarking that, with the rescaling Q(x) := Q{vqx ), we have: 


KnQ) = 

r o 


9 r 2 
z 'o 


|VQ| 2 + /(Q) 


dx + 


W 


2 r, 


0 J dB\ 


\Q s ~Q\ 2 dA. 


We will always work in this rescaled setting, and we assume from now on 
that ro = 1, and consider the domain 


fi := Hi = R 3 \ B, 


where B = B 1 is the ball of radius 1. The size of the particle is then encoded in 
the elastic constant L: L —> 00 represents the small particle limit, while L —> 0 
simulates very large particle radii. 

As mentioned in the introduction, an appropriate functional setting to es¬ 
tablish the existence of minimizers is the affine Hilbert space 



Note that the free energy functional (2) is not everywhere finite on the space 
Boo, since the potential term f(Q) > 0 may very well not be integrable in f2. 
However, since f(Qoo) = 0, we do know that 


inf T < 00 . 

«oo 

In fact, we show that the infimum is attained, and that the minimizer has a 
limit at infinity: 
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Proposition 3. Let L > 0 and W £ [0,+oo]. Then there exists Q £ P^ such 
that 

T(Q) = inf T. 

7~L oo 

Moreover, the far-field condition holds in the strong sense (5), and this is true 
for any solution of the Euler-Lagrange equation (10) in Poo- 

Remark 4. The case W = +oo will be understood as the strong anchoring 
case: it amounts to considering only maps Q £ Poo which satisfy the Dirichlet 
boundary condition (9) in the sense of traces. 

Proof of Proposition 3: Existence follows from the direct method of the calculus 
of variations. Thanks to Hardy’s inequality, any minimizing sequence ( Q n ) is 
bounded in Poo and admits (up to taking a subsequence) a weak limit Q £ Poo- 
We may also assume that the convergence Q n Q holds almost everywhere. 
Convexity and Fatou’s lemma allow us to conclude that P(Q) < lim inf P(Q n ) = 
inf T. 

The limit at infinity follows from estimates for solutions of the Euler-Lagrange 
equation (10). From Lemma 5 below we know that Q £ L°°(f2), and (10) read¬ 
ily implies that A Q £ Using standard elliptic estimates, we deduce 

that VQ € L°°(Ll), so that Q is uniformly continuous. Since on the other hand 
Sobolev inequality implies that \Q — Qoo\ belongs to L 6 (U), we conclude that 
\Q(x) — Qoo| converges to zero as \x\ goes to +oo. □ 

In the proof of Proposition 3 we used the following L°° bound for solutions 
of (10), related to the growth of the potential f(Q). 

Lemma 5. If Q £ Poo solves (10), it holds 

11*211 l°°{ n) — 9o> 

for some qo > 0 that depends only on a, b and c (but not on L and W). 

Proof. Let Q := Q — Qoo, so that Q £P solves 

LAQ = V/(Qoo + Q). (13) 

We may use as a test function in (13) any function T £ P with compact support 
in LI. Let us consider a test function 

U/ = VQ, for some V > 0 with compact support in H. 

Multiplying (13) by T we find 

L [ VQ • VT = - [ Wf(Qoo + Q) ■ Q + [ (14) 

Jn Jn Jon 

Clearly there exists q 0 = qo(a , b, c) >0 such that 

V/(<3oo + Q) ■ Q > 0 for \Q\ > qo, and \Q S ~ Qoo\ < Qo- 
Now let us take V of the form 

V = U<p, U = mm({\Q\ 2 -ql) + ,M), 0 <<p£C?(fi). (15) 
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Note that V is non-negative and supported inside the set {|Q| > < 70 }• 

Thanks to the choice of q 0 , both terms in the right-hand side of (14) are 
non-positive: for the first term this is clear, and the second term is zero in the 
case of strong anchoring (9) and non-positive in the case of weak anchoring (8) 
because \Q \ 2 — Q • (Q s — Qoo) > 0 for \Q\ > q 0 . Thus we obtain 

L I X7Q • V'k < 0, 

Jn 


i.e. 


j^\U\S/Q \ 2 + 
Next we take = ip R such that 


\m \ 2 


Pr(x) = 


for \x\ < R, 
for \x\ > 2 R, 


< - / UQ-VQ-Vp. 
Jn 


C 

and \Vip R (x)\ < —y, 


for some constant C > 0 independent of R. We obtain 


1 


PR ( U\VQ\ + g IVC/r ) < MC'||VQ|| L2( | x |> J e ) ||Q/r|| L 2 ( | x |> fl) . 


Since Q £ H, the right hand side converges to zero as R goes to +00 and it 
holds 

r (u\vq \ 2 


in 


V 


+ I |VC/| 2 ) = 0 


Recalling the definition (15) of U, we conclude that \Q\ < q 0 a.e., and therefore 

11*311 l°° — 9o + s*y / 2/3. □ 

Remark 6. It would be interesting to prove an explicit convergence rate for the 
far-field behavior (5). The small particle limit (cf. Section 3) suggests a bound 
of the form \Q(x) — Qoo| < C/ |x|, for some C = C(a, b , c, L ) > 0. An indication 
that this bound could indeed be true is given by the following proposition, which 
confirms that the minimizers Q approach the uniaxial set U* at the desired rate 
|a:| -1 . An analogous situation prevails in the case of solutions of the Ginzburg- 
Landau equations in the plane [27], which have an explicit rate of decay of the 
complex modulus « — 1, but provides much weaker asymptotic information in 
the behavior of the complex phase, the circle S 1 playing the role of 74 in that 
setting. 

Proposition 7. There exists a constant C = C(a,b 1 c) > 0 so that, for any 
solution Q £ RT° of (10) with finite energy T(Q) < 00 we have: 


C\/~L 2C 

dist((3(x),74) < and |VQ(a:)| < j^-r. 


(16) 


Proof. We follow the strategy of [27], which proved decay estimates for solutions 
to the Ginzburg-Landau equations in the plane. Let A := B 4 \ B\ and for any 
R > 1, A r := B 4R \ B r . Define Q R (y) := Q(Ry) for y £ A, with energy 

£r(Qr;A)= J A e R (Q)dy, with e R (Q)= ^ Q r\ 2 + ^ f (Q r) ■ 



By a change of variables in the integral, 

£ r (Q r -A) = ^J a ^|VQ | 2 + /(Q) dx = o(R~ 1 ), 

as R —> oo, by the finite energy assumption on Q. Thus, VQr —► 0 and 
R 2 I a /(Qr) dx —> 0 . Recalling that Q(x) —► Q^ as \x\ —> oo from Proposi¬ 
tion 3, we may assume Q R —> in H 1 (A). 

We now employ the convergence results for Landau-de Gennes as L —> 0, 
proven in [18], with our —> 0 replacing their L in this context. Although the 
convergence results in [18, § 4] are stated for global minimizers of the energy in 
bounded domains with Dirichlet condition, the proofs of the various convergence 
lemmas are based on the Monotonicity formula, and apply as well to solutions of 
the Euler Lagrange equations (10) with uniformly bounded energy, converging 
in H l -norm. In particular, we apply Lemma 7 of [18] to Q R in A := \ B\, 

which contains no singularities for sufficiently large R: since for every a £ A we 
have f Bl ( a j e B (Q B (y))dy < o( 1), by the Lemma there exists C 2 = 62 ( 0 , b, c) for 
which 

sup \1\VQr \ 2 + ^f(Q R )] <C 2 . 

Bi(a) Y l L J 

2 

Covering A by balls of radius 1, we obtain the same uniform bound 

lNQR(y)\ 2 + ^f(QR(y))<C2. (17) 

for all y £ A. 

Next, we claim that there exists a constant C 3 = Cs{a,b,c) such that for 
any Q G So, 

[dist(<3,74)] 2 < C 3 f(Q). (18) 

Indeed, since f(Q) > |Q | 4 as |Q| —> 00 and f(Q) = 0 if and only if Q G for 
any fixed 5 > 0 there exists C = C(S , a, b , c) such that the desired bound holds 
for all Q € S 0 with [dist(Q,W *)] 2 > 6 . 

Since LL, is smooth and compact, there exists 5 > 0 such that for all Q £ Sq 
with dist(Q,£A) < 6 there is a unique orthogonal projection Q* £ U * which 
minimizes the distance from Q to U*. Thus, Q = Q* + Q n , with Q n _L Tq„W*, 
and \Q n \ = dist(Q,W*). The function / is frame-invariant (it holds f(R~ 1 QR) = 
f(Q) for any rotation R), so without loss of generality we may assume that 
Q* = Qoo = S* (e z (g> e z — 1 1). Let us consider the orthonormal basis {Aj} of 

s 0 s iven b y 
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In this basis, TqJA * = spanjA^, A 5 }, and thus we may represent 

3 

Qn = Yu jAj . (19) 

i=1 

Expanding the potential in terms of u, 

f(Q ) = f(Q* + Q«) — f{Q*) = ^2a -I — ^ u\ H — («2 + ^ 3 ) + M U )> 

where h('u) = 0 (|u| 3 ) is a polynomial with terms of third and fourth degree in 
u. By the representation (19) of Q„, we thus have 

f(Q) = f{Q* + Qn) > Ci (u\ + v% + uf) = CilQn} 2 = C 4 [dist(Q,Z4)f , 

with constant C 4 = C 4 (S, a, b, c). The claim is thus established for all Q £ Sq. 

Finally, as in [27] we may return to the original scale x = Ry, to obtain the 
desired bounds (16). Indeed, by (17) and |x| = i?|y| < 3 R, 

i|VQ(x)| 2 + I/(Q(*)) = R- 2 Q|VQ fl (2/)| 2 + ^.f(Q R (y))j < g < p, 

so the conclusion follows from the above and (18). □ 

Remark 8 . From the decay estimate (16) it is straightforward (see Proposi¬ 
tion 7 of [18]) to show that the eigenvalues of Q(x) tend to those of the uniaxial 
tensors in W*, and at the same asymptotic rate \[L/\x\ as |x| —> oo. 

3 Small particle: the Saturn ring 

This section is dedicated to proving Theorem 1, and then studying the limiting 
configuration Q 0 , whose expression we recall here: 

Qo = s* » W 4 (e r 0e r - \l ) + s*(l - “) \ e z®e-z - ) • ( 20 ) 

3 + w r A \ 3 J 1 + wr\ 3/ 

Proof of Theorem 1. Recall that we assume r o = 1, so that we are in fact taking 
the limit (1 /L,W/L) —» (0,iu). It is straightforward to check that the map Qo 
( 20 ) belongs to Hoc and solves 


AQ o = 0 in Q, 

1 d Qo o co oo 
-= Q s - Qo on dil. 


( 21 ) 


In what follows we emphasize the dependence of Qo on the parameter w > 0 by 
writing Qo = Qo,iu- Consider the map Q := Q — Qo,w/l, which solves 


j AQ = I-V/(Q) in 0, 

1 TTr TT + Q = 0 on 
v W ov 
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Next we apply an analog of the interpolated estimate of [3, Lemma A.2] for the 
oblique derivative problem (boundary conditions du/dv+au = 0) in fl: it holds 

fjVQll^ <C||AO|| ioa ||Q|| Xoo . (23) 

This estimate can be proven exactly as in [3]: apply [3, Lemma A.l] which does 
not need to be adapted, and then follow the proof of [3, Lemma A.2] using 
elliptic L p estimates for the oblique derivative problem near dfl (see e.g. fl, 
§ 15] or [6]) instead of the homogeneous Dirichlet problem. The constant C > 0 
depends only on <5 > 0 such that w > 5. 

Since we can infer from Lemma 5 that ||V/(Q)|| i0 o < C (where C > 0 
depends only on a, b and c), the estimate (23) and the equation (22) imply that 

IIVQlLoc < (24) 


from which we deduce that for any compact K C f1, it holds 

diam(A') 




Vl 


with the constant C > 0 depending on a, b, c and 6 > 0 such that w > 5. On 
the other hand it can be easily checked that, at any order k € N, 

||Qo,W/L — Qo,u)|| C fc < Ck |l/w — L/W | , 


so that Q — Qo = Q + Qo,w/l ~ Qo,w converges to zero locally uniformly in fl, 
as (1 /L,W/L) (0, w). ' ' □ 

The main interest of Theorem 1 lies in the explicit expression for the limit 
Qo- Next we investigate its most important features. We start by interpreting 
the Saturn ring defect as a locus of uniaxiality. More precisely, let U w be the 
uniaxial locus of Qo = Qq w (20) away from the 2 -axis (on which e r = ±e z and 
Qo is trivially uniaxial): 


U w := {x € fl \ Re z : Qo(x) is uniaxial} . 


(25) 


Then we have: 

Proposition 9. For w > \/3 it holds 


U w = {(aq, x 2 , 0): x\ + x\ = r^} , 


where r w is the unique solution r > 1 of 


r 


3 


w 9 w 

- r - 

1 + w 3 + w 


(26) 


(27) 


The function w <—> r w increases from r^ = 1 to a finite value r^ ~ 1.47, as w 
increases from -\/3 to oo. 

For w < v/3, U w is empty. 


11 



Proof. Using spherical coordinates (r, 0,ip), the map Qo is of the form 
Qo = a(r)(e r <g> e r - 1/3) + fi(r){e z <g> e z - 1/3), (r > 1), 
where the functions a(r),/3(r) > 0 are given by 


V» = ®*(! TU~ ~)' 

w r 6 1 + w r 


Using the fact that e r = cos^j e z + sin<p e p , and defining 


a := a + /3, v := a/3, 


(28) 


elementary computations show that the characteristic polynomial of Qo is 



so that the above square root is well defined. The eigenvalues A® > A° > A° of 
Qo are given by 



We are looking for the points where Qo is uniaxial, i.e. either A® = A§ or 
A° = A°. For 0 < ip < 7r, it holds 

A? = A 2 p = ^ and a = /3, 

A° = A° 4==> a = 0 or /3 = 0. 

Since a,/3 > 0, only the first case (A? = A°) can occur. Given the expressions 
(28) of a and /3, we deduce that 

u w = {(x!,x 2 , 0 ): 1 < x\ + x\ = r 2 , p{w,r) = 0} , 

1 r \ 3 w 2 w 

where p[w,r):=r - r -. 

1 + w 3 + w 

It is straightforward to check that r 1 —> p(w,r) is increasing on (l,oo) and 
therefore p(w,-) = 0 has at most one solution in this interval. Since on the 
other hand 

(1 + w )(3 + w) ■ p(w, 1) = 3 — w 2 , 
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there is a solution r w > 1 if and only if w > \/3. The function w K > r w is easily 
seen to be smooth, and its derivative has the same sign as 


—d w p(w, r) 


(1 + w ) 2 


3 

(3 + w) 2 


> 0 , 


so that r w is an increasing function of w. As w increases to oo, r w increases to 
Too > 1 such that — 1 = 0. □ 

Remark 10. Since Q o is biaxial except along the Saturn ring locus U w , there is 
no conventional nematic director field attached to it. However, it is reasonable 
to distinguish a principal direction as an approximate or mean “director ’’field no, 
the unit eigenvector associated to the largest eigenvalue A?. This is well-defined 
up to a sign, at every point where is a simple eigenvalue: that is, everywhere 
except at U w , where it jumps discontinuously as the eigenvalue branches cross. 
Then one may compute, for 0 < ip < 7t/2, 


n 0 = 



T := 


a(l — 2 sin 2 ip) + p 


a 2 + p 2 + 2a/3(l — 2 sin 2 ip) 


with a, P as in (28). For it > <p > it /2 the director field is obtained by reflecting 
with respect to the horizontal plane : n 0 (<p) = y^l — p )/2 e p — \/(l + p)/2 e z . 
Note that this way no is not continuous in H \ U w , but it is continuous as an 
MP 2 -valued map : the map no ® no is continuous in Q \ U w . 


Corollary 11. Let w > \/3. There exists <?o > 0 such that, if 

r 1 L 1 

^ -J + 777 -< 

L W w 

then any solution Q £ TLoo of (10)-(8) admits a uniaxial ring near U w . More 
precisely, Q(x ) is uniaxial for all x belonging to 

U = {{pu(0)cost),p u (6)sin0,z u (0)): 9 e K} , 

where the functions p u (9) and z u (9) satisfy 

I Pu(9) - r w \ + \z u {6)\ < e(5) 0, 


as S —> 0. 


Proof. Biaxiality can be quantified through the biaxiality parameter [15] 

,(MQ 3 )) 2 


P{Q) = 1 - 6 - 


\Q\ 


which is such that : Q is uniaxial if and only if p(Q) = 0. Let us fix R > r w . 
From Theorem 1 and Proposition 9 we infer that there exists (j 0 such that for 
6 < S 0 it holds 


P(Q) >0 in A := {x £ Br : dist(x, U w U Ke 2 ) > e(<5)} , 
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(a) w = oo. 


(b) w = 3. 




Figure 1: The director field no (integral curves). 


for some e( 8 ) —> 0 as 5 —> 0. Let us write xo, the characteristic polynomial of 
Qo, as 

Xo = (X - A?)P 0 , 

where Po = (AT — Aj^X — A 3 ). By continuity of the roots of a polynomial, it is 
clear that the characteristic polynomial % of Q satisfies 

X = (X - Ai)P, |Ai - A?| + |P - P 0 | < ci(5) in A, 

for some ci(5) —> 0 as S —> 0. The eigenvalue Ai and the coefficients of P depend 
continuously on x £ A. Note that Po(Qo)no ^ 0 in A, and define 

“ := iwkii p(Q) ’‘" 

Then Qu = Aiit, and |« — no| < 02 ( 8 ) —> 0 in A, so that 11 / 0 and we may 
define n = u/ |u|. It holds Qn = Ai?r and 

| n{x) - n 0 (x )| < c 3 (S) -> 0. 


14 







Moreover the map n ® n is continuous in A. 

Now fix 0 G K. and denote by He the half plane corresponding to the az¬ 
imuthal angle 6 


He = {{p cos 6 , p sin 9, z) : p > 0, ztR}, 

and by De the disc in He of radius e(<5) and of center at the point where the 
ring U w intersects Hg: 

Dg = {(p cos 9, p sin 9, z): (p — r w ) 2 + z 2 < e(c>) 2 } . 

We claim that for low enough 6 there exists x G De such that Q(x) is uniaxial 
(which obviously proves Corollary 11). 

To prove the claim, note that when restricted to He, the director field no 
may be viewed as an S 1 -valued map since it takes values in S 2 ft He- Then the 
restriction of no ® no to dDg is topologically non trivial : it corresponds to a 
non trival class of 7 t 1 (MP 1 ), as can be seen by explicitly computing its degree. 
Since for 5 low enough n{x) is arbitrarily close to no (a;), this implies that the 
map n®n, which is continuous in AtlHg, admits no continuous extension inside 
De- Therefore Q can not be biaxial everywhere in De : if it were the case, Q 
would have only simple eigenvalues and admit a differentiable eigenframe [ 21 ]. 
In particular there would be a differentiable vector field n defined in De, such 
that h ® h extends n® n. □ 


4 Large particle: the dipole structure 


As mentioned in the introduction, the core of Theorem 2 is the fact that the 
axially symmetric harmonic map obtained in the limit of a large particle has 
exactly one singularity. In this section, we prove this result (Theorem 12 below) 
and then complete the proof of Theorem 2. 

We define the axially symmetric S 2 -valued maps to be exactly the maps 
n G Hl oc (Q; S 2 ) which can be written in cylindrical coordinates (p,9,z) as 

n = smtjj(p, z) e p + cosip(p, z) e z , (29) 

for some real-valued function i)) G Hi (Q. cy {) defined in the domain 

^cyi := {{p,z) G M+ x R: p 2 + z 2 > 1} . (30) 

We consider here strong anchoring conditions given by 


n = e r for |x| = 1 , 

and the far-field conditions in integral form 

f (ni ) 2 + (n 2 ) 2 


dx. 


(31) 


(32) 


in x 


As in Section 2, the existence of an axially symmetric S 2 valued map n mini¬ 
mizing the Dirichlet functional 


E(n) = [ |V?r | 2 , 
Jn 
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under the conditions (31)-(32) follows from the direct method of the calculus of 
variations and Hardy’s inequality. Such a map is analytic away from a discrete 
set of singularities on the 2 -axis [9]. Here we prove: 

Theorem 12 . Let n £ H[ oc (Q; S 2 ) be a minimizer of the Dirichlet functional 
E among all axially symmetric S 2 -valued maps satisfying (31)-(32). Then n is 
analytic away from exactly one point defect: there exists |zo| > 1 such that n is 
analytic in \ {( 0 , 0 ,2 o)}. 

Proof of Theorem 12. Preliminaries. Let if £ Hl oc (^cyi) be the function as¬ 
sociated to n through (29). Then if minimizes the energy 


E(if) = 



\d p if \ 2 + \d z if\ 2 + -4 sin 2 if 
P 2 


pdpdz, 


among all functions if £ Hl oc (fl cy i) such that the corresponding n satisfies (31)- 
(32). Next we express these boundary conditions in terms of if. 

The strong anchoring condition (31) is more conveniently expressed using 
spherical coordinates (r,9,ip): 


if = ip for r = 1. 


(33) 


It holds in the sense of traces, which makes sense as soon as E(if) < oo. (In 
fact we should have written if = ip mod 27 t, but since any Z-valued function of 
regularity ff 1 / 2 is constant [4] we may reduce to the above.) 

The far-held condition (32) becomes 



sin 2 if 
p 2 + z 2 


pdpdz < oo. 


(34) 


Hence the class of admissible functions consists exactly of the if £ Hf )c (£l cy [) 
satisfying E(if) < oo and (33)-(34). 

The Euler-Lagrange equation satisfied by if is 

d 2 z if + d 2 p if + -d p if = sin(2-0) in f l cyi . (35) 

P ^P 

Note that, by elliptic regularity, any solution of (35) is real-analytic away from 
the 2 -axis {p = 0}. Also note that, since replacing if by max((), 0) or min(^, tt) 
does not change the boundary conditions and decreases the energy, it holds 


0 < if < TT. 


The rest of the proof is divided into 3 steps. 

Step 1. We claim that the open subsets of fl cy i n {p > 0}, 

X + = {if > 7 t/ 2 } and = {if < 7 t/ 2 }, 


are connected. 

We split the half-circle dfl cy i D {p > 0} into two arcs: 

A± = {r = 1, ip = tt/2 ±t: t £ (0,ir/2)} . 
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The boundary conditions ensure that A± C X±. We denote by u>± the connected 
component of X± containing A±. 

Let us show first that X + = ui+. Consider the function 




ip in w+, 

nrin('0,7r — ip) in f l cy i \ u>+. 


Then it can be checked that ip G Hl oc (Cl cy i). Moreover, E(ip) = E{ip) and ip 
clearly satisfies the boundary conditions (33)-(34). Therefore ip minimizes E, 
and is analytic away from the 2 -axis. In particular, min(^>, 7 t — ip) is analytic in 

cyl \ • 

On the other hand, since X _ is open and non-void, there is an open subset 
of f l cy i \ w+ in which ip < tt/2. In that open subset, the two analytic functions 
ip and min (ip, n — ip) coincide, so they must coincide in the whole Q cy i \ w+. We 
deduce that 

ip < tt — ip i.e. ip < 7t/ 2 in Q cy i \ u;+, 

and therefore X + = is connected. 

To show that X_ is connected, consider 


ip 


iP 

ma x(ip, tt — ip) 


in w_, 

in f l C yi \ w_. 


As above, E(ip) = E(ip) and we conclude that AT = w_ is connected. 

Step 2. There is at most one singularity. 

We know [9] that n is analytic in fl away from a set of isolated points 

Zc{p = 0, \z\ >1}. 

In particular ip is continuous in Q cy i \ Z , and since 

[ sin _4p dpdz < < 0Q) 

J fl C yl P 

and 0 < ip < n, it follows that 

ip G {0, tt} on ( dn cy i n {p = 0}) \ z. 

At every point (0, zq) G Z, ip must be discontinuous, because otherwise n would 
be continuous around that point (and then real analytic). Therefore it must 
hold 


either ip = 


0 in (zo - S,z 0 ), 
tt in (z 0 , z 0 + S), 


or ip = 


tt in (20 — S, Zq), 
0 in ( 20 , 2:0 + 6), 


for some S > 0. 

Let us argue by contradiction and assume that there exist two distinct points 


(0, 21 ), (0, 2 2 ) G Z, 21 < 22 , [ 21 , 22 ] <1Z = { 21 , 22 }. 

There are three cases: either 21 <22 < —1, or 21 < — 1 < 1 < 22 , or 1 <21 < 22 . 
Note that the boundary conditions (together with the boundary regularity of n 
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[9]) ensure that ip = n in (—1 — 5 , —1) and ip = 0 in (1,1 + 5) for some S > 0. 
In all three cases, it is easy to see that there must exist four distinct points 

( 0 , z'j) £ Z, z[ < z 2 < z 3 < z 4 , 

such that 

either ip(z[) = ip(z' 3 ) = 0, ip(z 2 ) = ip(z' 4 ) = n, 
or ip(z[) = ip(z 3 ) = 7r, ip(z 2 ) = ip(z' 4 ) = 0. 


We may assume that we are in the first case (the second case can be dealt with 
similarly). Then by continuity there exists 8 > 0 such that 

[B(( 0 , 4), 6) n U cyl ] u [B((o,4), 6) n U cyl ] c a_, 

[B((0,4). 4 n n cyi ] u [B((0,4). 6) n Ti cyl ] c x+. 


Since the sets X± are path-connected we deduce from the above that there is a 
continuous path from z[ to z 3 inside X_, and a continuous path 7 + from z ' 2 
to 4 inside X+. Since 4 < 4 < 4 < 4 > these paths must intersect, but then 
there intersection would belong to DX + = 0. This contradiction shows that 
Z contains at most one point. 

Step 3. There is at least one singularity. 

Assume that Z is empty. Then n is continuous in f 1, and therefore 


deg(n, dB r ) 


dcpip(r, ip) sin ip(r, ip) dip, 


is independent of r > 1. We deduce that 


1 = deg(n, dB\) = deg(?r, dB r ) 

1 
4 

which implies E(ip) > 4 dr = 00 , a contradiction. □ 

Before turning to the proof of Theorem 2, we define rigorously the axially 
symmetric Q-tensor maps. They are the maps Q £ Boo which satisfy the two 
following natural symmetry constraints: 

• the map Q is invariant by rotation around the vertical axis: 

Q(Rx) = t RQ(x)R for all rotations R of axis e z , (36) 


d v ip{r,ip ) 2 sin <pdip, 


1 

< - 
“ 4 


sin 2 ip(r,ip) . 


sin 2 ip 


sin ipdip ■ 


• the vector eg is everywhere an eigenvector of Q: 

Q(x)eg ■ e p = Q(x)eg ■ e z = 0. (37) 

These constraints are natural, in the sense that a minimizer of the free energy 
( 2 ) under those restrictions is still a solution of the complete (unconstrained) 
Euler-Lagrange system (10) - as can be easily checked. We denote this class of 
symmetric maps by B s ^ m C Boo- 
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Proof of Theorem 2. Recall that we assume ?’o = 1 and we are therefore consid¬ 
ering the limit 1/L —> oo. Since axial symmetry (36)-(37) is clearly preserved by 
pointwise convergence, we may proceed exactly as in [3, Proposition 1] to obtain 
a subsequence (Qk) converging in %oo to an axially symmetric W*-valued map 
Q* minimizing the Dirichlet energy j n |VQ| . The estimates in [18, 20] show 
that the convergence is in fact locally uniform in fl, away from the singular set 
of Q*. 

Because fl is simply connected, ZT-valued Hl oc maps can be lifted to S 2 - 
valued Hj oc maps [2]: for any 7/*-valued Q £ Hoo, there exists n £ H/ oc (fl; S 2 ) 
such that 

Q Qn •— S* (^ ^ ri I /3) ■ 

Therefore, the limiting map Q* can be written as Q* = Q n , and the map n 
minimizes the Dirichlet energy in the class 

'Hl ym := {n £ Hj oc (fh § 2 ): Q n £ R s ^ m with strong anchoring (9)} . 

To conclude the proof, it remains to show that the class TLt vm actually 
corresponds to the class considered in Theorem 12. 

Using the fact that |n| 2 = 1, we calculate 

s* 2 |Qn ~ Qoo| = 2(n 2 + n\), 

so that Q n satisfying (11) translates into n satisfying (32). 

The strong anchoring condition (9) for Q n is equivalent to 

n \dB = T e r, 

for some {±l}-valued function r, which must be of regularity if 1 / 2 and therefore 
constant [4[. Therefore, up to multiplying the map n by a sign, the strong 
anchoring condition becomes (31). 

The fact that Q n admits eg as an eigenvector (37) is equivalent to 

(n ■ eg)(n ■ e p ) = (n ■ e e )(n ■ e z ) = 0 n • e e G {0, ±1}. 

Since the function n ■ eg is flj 1 , it must therefore be constant. The boundary 
conditions prevent it to be equal to ±1, so that n ■ eg = 0. 

The invariance by rotation (36) for Q n is equivalent to n(Rx) = ±Rn(x) for 
any rotation R of axis e z . The sign ±1 may depend on x and on R, but the 
Hloc regularity implies that it does not depend on x. Therefore it holds, using 
cylindrical coordinates, 

n(p, 9 , z) = t(9) Rg n(p, 0, z), t{9) = ±1, 

where Rg £ 50(3) is the rotation of axis e z and angle 6 . The function r is 
easily seen to belong to T7 1 / 2 (K/27rZ): we conclude that r = 1. 

Therefore the axial symmetry (36)-(37) of Q n is equivalent to 

n(p,9,z) = Rgn(p,0,z), and n ■ eg = 0, 

which implies that 

n(p,6,z) = u 1 (p,z)e p +u 2 (p,z)e z , u £ Hl oc (fl cy i; S 1 ). 

Since fl cy i is simply connected, u can be lifted to a real-valued function ip £ 
HlodilcyO: u\ = sin^>, it2 = cos ip. Therefore n is of the form (29), and Hi vm 
corresponds indeed (up to a sign) to the class of axially symmetric S 2 -valued 
maps satisfying (31)-(32). □ 
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